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ABSTRACT 


A simple system was designed to collect wide-angle reflection 
records in order to investigate the interval sound speeds of the sedi- 
ment layers in the ocean. The system consisted of a frequency 
modulated receiver, a cut-to-channel Yagi antenna, and a sonobuoy, 
used in conjunction with a precision sonic profile recorder and a 
triggered sound source. A computer routine for reducing the data was 
obtained and modified for compatibility with the system and use on the 
IBM 360/67 computer. The system was designed to be both inexpensive 


and simple to use without any loss of accuracy. 
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Pe INTRODUCTION 


The objective of this study was to develop a simple yet in- 
expensive system for collecting wide-angle reflection data. This 
system would be used in conjunction with normal-reflection profiling 
equipment normally carried on board oceanographic survey vessels. 
The study involved testing an inexpensive radio receiver to be used to 
receive a signal transmitted by a sonobuoy with enough clarity that 
the signal could be recorded on a precision seismic graphic recorder. 
The study also involved testing and launching five sonobuoys and in- 
vestigating modifications which can be made to the sonobuoys to 
enhance the success of the system. 

A further objective of the study was to provide a method for re- 
ducing the data collected by the system and obtaining the mean sound 
speed existing between successive reflectors (called the interval 
sound speed of the layer defined by an upper and lower reflector) and 
the thicknesses of the corresponding sub-bottom layers. This required 
modifying an existing program to read the data and compute the results 
on the IBM 360/67 digital computer using FORTRAN IV language. 

This paper presents a short discussion of the background and 
purpose for developing a system for collecting wide-angle reflection 
data and obtaining a program to reduce the data and determine the in- 
terval sound speeds. Atheory of wide-angle reflections is offered, 


together with a mathematical model used to develop the equations and 





suggest a method of solution. The computer program is discussed, 
including the method of digitizing data from a wide-angle reflection 
record. The radio receiver, antenna, and radio sonobuoy used to 
develop the system for use in conjunction with normal-reflection 
equipment is discussed. In the last section, conclusions and sug- 


gestions for further work on the subject are offered. 


10 





II. BACKGROUND 


The geological structure and composition of the sub-bottom 
layers in and around Monterey Bay are presently known to only a 
limited extent. Sufficient depth recordings have been made to enable 
drawing a detailed relief map of the sea floor; and shallow cores and 
grab samples have been analyzed for determination of the composition 
of the sediments in the first few meters. But knowledge of the deeper 
layers has been restricted to only subjective analysis involving the 
extrapolation of the immediate terrestrial structure from wells and out- 
crops and correlating the data with normal-reflection seismic profiling 
data. Unfortunately, normal-reflection profiles measure only the 
travel times and relative intensity of the vertical reflections. Without 
knowledge of the sound speed-depth relationship, only a guess can be 
offered as to the actual sub-bottom structure. A knowledge of the 
sound speed-depth variation is essential to a quantitative evaluation 
of the sediment layers and to a geological interpretation relating sound 
speed to type of material, age, and consolidation (Le Pichon, Ewing, 
and Houtz, 1968). 

The system developed for collecting wide-angle reflection data 
along with the computer program to reduce the data can provide in- 
formation on the sound speed-depth relationships and the layer thick- 


NEessSes . 


et 





Pie die Oo-ANGLE REFLECTION METHOD 


The characteristic feature of the wide-angle reflection method is 
the measurement of the travel times of longitudinal waves which have 
been reflected at boundaries separating media of different acoustic 
iaredances (the product of sound speed and density). From measure- 
ments of reflection times it is usually possible to determine the depths 
and dips of the reflection horizons and the speed of the seismic wave 


(JIakosky, 1949). 


A. Peet ORY 

The system was developed based on the theory that a sound pro- 
duced at source S (Fig. 1) will reach the hydrophone of the sonobuoy B 
by three different ray paths as shown; ] the direct ray which travels 
just below the surface of the water, 2 the sound reflected from the 
water-sediment interface, and 3 the sound waves reflected from the 
various sediment layers that may exist below the ocean floor. If the 
distance separating the sound source and the hydrophones is steadily 
increased while the source emits sharp sound pulses at regular intervals, 
the time difference between the arrivals of the direct and reflected 
sound waves will change depending on the sound speed-depth variation 
within the sediments. A record of time differences is used to compute 
the thicknesses and interval sound speeds of the various layers and 


‘obtain the sound speed-depth relationship. 
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re THE MATHEMATICAL MODEL 
le orming the T/X Plot 

For the purpose of forming an elementary concept of the 
reflection method and to illustrate the formation of the wide-angle re- 
flection records, it is useful to first consider a simple case ofa 
single layer of homogeneous isotropic material of uniform thickness H. 
It follows from this case that the ray paths will follow straight lines. 

A further assumption may be made that the lower surface reflects the 
sound wave in accordance with the physical laws of optics so that the 
angle of reflection is equal to the angle of incidence. Also it will be 
assumed that the path taken by sound passing from One media to another 
having a different sound speed is refracted in accordance with Snell's 
Law. 

Figure 2 illustrates a single layer model with a water 
medium of depth H and average vertical sound speed V. The theoretical 
travel time on the record for a vertically-propagating sound wave when 
the sound source and sonobuoy are at the same point S in the ocean is 
referred to as TO: Thus To defines the origin of the record at which 
point the direct travel time is zero. The source sound wave will 
propagate downward through the water to point A and then be reflected 


back to the sonobuoy in time 
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Figure 2. 


reflected from a horizontal surface 
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If an anology is made to the laws of geometric optics, the 
ray paths may be conveniently handled by utilizing the concept of 
images. The travel time to the theoretical image point 0’ is equal to 
the time for the reflected ray to return to the sonobuoy. Therefore, 
when a mathematical model is constructed, the reflected times can be 
considered as the time required for the sound to travel from any shot 
point 5. to the image point 0' at a sound speed of V. 


At shot point S,, a sound wave travels to point B, where 


] l 


it is reflected upward to the sonobuoy. This wave is referred to as 
the reflected ray from the first reflector. The direct ray travels hori- 
zontally from the source to the sonobuoy at S. Again using the concept 


of images, the reflected ray path S8.B.S is equal to the assumed path 





: i 
5B, 0 (path length R,) and the travel time becomes 
Be 
= ? 
fe v7 ; (2) 


where qT is the travel time of a reflected ray. From the geometry of 


the triangle SBA 


2 2H 
ro (3) 
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Equation (2) then becomes 


- 2H 
i Vcos#, (4) 
where B, is the angle of incidence measured from the vertical. 


The direct ray arrives at the sonobuoy at time 


where D, is the travel time of a direct ray, X 


j is the separation 


] 
between the shot and the sonobuoy, and V is the sound speed. 
In general, the separation X can be related to the depth H 


and the length of the reflected ray path by use of the Pythagorean 


Theorem; i.e., 


where R is the distance from the shot point to the image point 0’. 

As successive shots are fired while steadily increasing 
the separation, the record of the direct travel times results ina straight 
lnievometne 17% plot shown im the upper half of Figure 2. From Equation 
fomiutccanibe seen that the slope of this line is the reciprocal of the 
water sound speed. The record of the arrival of the reflected rays form 
a hyperbola asymptotic to oe direct travel time line at large values 


ly 





of X. Equation (6) can be combined with the relation 


where T is the travel time of the reflected ray and V is the sound speed 


to obtain the equation of a hyperbola 


(7) 


For a bottom with a slope shown in Figure 3, it can be 
shown that an additional term is introduced into the above equation. 
The equation for the reflected travel time for such a model is 

] 2 7 


T= ee a + 4HX sin 8) 


e (8) 


where 6 is the angle of slope of the bottom. The sign of the additional 
term is negative if the bottom dips in the direction opposite to that 
shown in Figure 3. 

If a mean sound speed between the sea surface and any 
sub-bottom reflector is considered, it will be realized that the record 
of reflection times from that reflector will also form a hyperbolic trace 
on the T/X graph in accordance with Equation (7), with V equal to the 
mean sound speed between the surface and the reflector and H equal 
to the depth to the reflector. The mean sound speed V can be calculated, 


but it provides little information for determining the sound speed-depth 
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relationship. The sound speed that is required for determining this 
relationship, however, is the interval sound speed of the layer defined 
by an upper and lower reflector. 

Zi. DeTeanination of Interval Sound Speed 

The determination of the interval sound speeds by reflection 
measurements is an old geophysical problem (Clay and Rona, 1965) and 
techniques for reducing the data go back a number of years (Green, 1938). 
The method of determining the interval sound speeds between successive 
reflecting layers from wide-angle reflection data is presently quite well 
known and has been described extensively in the literature (Dix, 1955). 
Its use at sea has been discussed by Houtz and Ewing (1963), Knox 
(1965), Clay and Rona (1965), and Le Pichon, Ewing, and Houtz (1968), 
among others. The method used in this paper follows the technique 
introduced by Le Pichon, Ewing, and Houtz (1968) and is the method 
used in the computer program of Appendix A. 

The water layer is assumed to have a known sound speed- 
depth profile, with the horizontal sound speed at the surface Vi dif- 
ferent from that of the sounding speed (the sound speed integrated 
over the depth H of the water layer). 

It is further assumed that the dip of the bottom is known. 
The dip of any sub-bottom layer will be considered here to be the angle 


between the upper and lower reflectors defining the layer. Dip is con- 


sidered positive if the distance between the two reflectors increases 


20 





in the direction of separation between source and receiver. The dip 
of the water layer illustrated in Figure 3, for example, is positive. 
Assuming such a dip and using a mean water sound speed 
equal to the sounding speed oe Equation (8) can be rewritten as 
T° = se x? + aH + 4HXsinB), (9) 
Mi 


and Equation (1) becomes 


eee age, ee (10) 


Since the direct travel times are a function of the hori- 


zontal sound speed V, , Equation (5) becomes 


h 


eee, Ge) 


Substituting Equation (10) and (11) into Equation (9) and 
rearranging terms, an equation relating the reflected travel time to the 


direct travel time is 








) Dsin @. (12) 


If Vv, is replaced by the mean vertical sound speed of the material con- 
tained between any two reflectors, called the interval sound speed 


ea , Equation (12) will relate the reflected travel time within the 
layer to the direct travel time at the surface corresponding to the 


Zt 





reflected time. The subscript i+] indicates the sound speed for the 
i-th sub-bottom layer. The wide-angle reflection records, however, 
relate total reflected travel time to total direct travel time. 

The technique used to solve for the interval sound speed 
involved reducing the travel times for effects of the upper layers in 
order that the reduced values have the form of Equation (12). An 
iterative method is then used to reduce the effect of dip by eliminating 


mafesecunvature term 8 Of Equation (12): i.e. , 





Finally, a least-squares line is fitted to the reduced 





2 Z 
(T= T)/D° data and the coefficient 
on 2 
( vy ) 
ictede 


is determined. From this, the interval sound speed is readily 
determined. 
aie First Layer 
The water between the surface of the ocean and the 
sea floor is considered here to be the first layer. It is assumed that 


the dip of the sea floor 9, is known and the sounding speed V4 can be 


1 


determined with good accuracy. 
Knowing aa and 8, and knowing the reflected travel 


times as a function of the direct travel times from the T/X data, the 


Ze 





solution of the first layer involves computing the thickness of the water 
layer Hy and the horizontal sound speed at the surface va? 
If data are obtained close to Toe a fourth-order 


least-squares polynomial is fitted to the T/X data, 


2 5 “ 
=a + +f + + 
It ay a,D aD a,D a,D ; (14) 


where the a's are coefficients of the polynomial. At reflected travel 
time, ee the direct travel time D is zero; the reflected travel time, 
also being the minimum travel time for the reflected wave, is equal 


to the coefficient ani inex, 
f= ate: (15) 


The fourth-order curve is then extrapolated, obtaining the minimum 
reflection time ie 

Care must be used with this method because the 
least-squares fitting technique does not constrain the fourth-order 
curve beyond the limits of the data. A more reliable but less accurate 

ae. ; 
method is to fit the T /X data to a linear least-squares line of the form 
2 2 


£f =bb.F b-D . (16) 
O 2 


. Z 
Again it can be shown that the coefficient dD, is equal to ue » hase 
Ty is obtained by extrapolating a straight line instead of a curve as 


was done in the case of the fourth-order polynomial. 
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The computer program used to reduce the TX data 
decides which of the above two methods will be employed to solve for 
Le Three-tenths of a second direct travel time was considered the 
maximum allowable distance from which the curve could be extrapolated 


using the fourth-order polynomial. 


Having computed Ty and knowing V., the depth of the 


1 


leer a- is computed Using Equation (1). 


1 
The curvature term in Equation (12) is small compared 
to the other two terms and, therefore, little error is introduced by 


assuming V/Vy equal to unity. This assumption is used to define a 


correction term to be applied to Equation (12); i.e., 


bat? +27 (1) sin @ (17) 
O O 


The reduced travel times T are then computed for each data point using 
iG 


the relation 


ee ae ; (18) 
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where T. is a reduced travel time with the effects of both the dip and Ty 
eliminated. The horizontal water sound speed Mh fssthen found by hitting 


Z 2 
rT and D to a least-squares first-order line of the form 


T =c.D, (20) 


emeesolving for the coefficient Cy: Knowing Vy , Equation (19) can be 


solved for Vii ey 


The separation X is then computed using Equation (5), 


a (22) 


De Second Layer 

The second layer is, in fact, the first layer of sedi- 
ments for which an interval sound speed is desired. The record of 
reflected travel times from the lower interface of the second layer qT, 
forms a hyperbola similar to that of the first layer. In the first layer 
the interval sound speed yy was known and the separation X was com- 
puted for each point. In the second layer, the interval sound speed V, 
is unknown, but the reduced direct travel time can be deduced by reducing 


the effects of the first layer on the travel times of the second layer. 


25 





The thickness of the second layer H, is computed 


2 


and dip oe) of 
a a 


using assumed values for the interval sound speed Vy 


the layer. These values are referred to as the trial solution for the 
second layer. The trial solution is based on the apparent dips read on 
the normal reflection profiles, converted to true dips using trial sound 
speeds. The exact value of these trial sound speeds are not critical, 
so long as they are consistent with the dips entered. 

In Figure 4, an assumed ray path is drawn as the 
dashed line SBC. This ray is perpendicular to the second reflector; i.e., 


to the lower interface of the second layer. From Snell's Law 


V 
V, sin 8.) : (23) 





sing= ( 


The distance the ray travels in the first layer SB is 


therefore 


= Ay 
Bese 24 
B cos g ‘ ( ) 


and the time Ty in the first layer 


(25) 


where MG and Hy are both known. 


200 





SEA SURFACE . 





Figure 4. To reflection from the second layer 


2. 





To obtain the reduced minimum travel time ET for the 
2 


r 
second layer alone, the minimum measured time to the second reflector 


Ty (obtained fitting the data to a least-squares polynomial) must be 
2 


reduced by the time the sound traveled in the first layer Te 


The length BC is obtained using the trial sound speed 


BC “We er ; (27) 
a 2. 


The length AB is the distance the direct ray will travel 
corresponding to the distance SB the reflected ray travels in the first 


layer. AB is computed from the expression 


AB =H, tan Z, (28) 


and the depth of the second layer H, below the shot point S at reflected 


2 


travel time T, is 


ca 
— 


= , 29 
H, B + AB sin 8, (29) 


Hi, is a first approximation for the thickness of the second layer. 


28 





Since V, was computed for the first layer, Equation (14) 





h 
can be written as 
=a x + eo ae ae x4 (30 
O 1 2 aS 4 ’ 
where 
X=DV., 
h 
and 
a, 
_ i 
oe Vv, 


' Differentiating with respect to X, the horizontal 


separation, Equation (30) becomes 


oF 


T 2 3 5 
ry = dy 2dx 3d 4X + 4d 1 P (3 1) 


It can be shown (Clay and Rona, 1965) that the angle of emergence 


at the sea surface §+6, is related to the derivative by 


sin (g + 8,) 


mre (", (32) 
Mi 


S| 


The data from the T/X graph is fitted to a fourth-order 


least-square polynomial and the coefficients of Equation (30) determined. 


Zo 





The angle (@+@) is computed for each data point by combining Equations 


(o) and (32) to give 


2 3 
t + = + 
sin (g 8,) V, (d, 2d,X + 3d.4X + 4d x ve (33) 


Knowing the dip of the first layer 8., @ is easily obtained. 


1 f 

From %, and using the sounding speed We , the travel 
time of the second reflected ray within the first layer is computed. This 
is then subtracted from the total reflection time to obtain the reduced 


metlectea trave! time T,, Pine time the retlected sound from the second 
i 


reflector travels within the second layer. Similarly, a reduced direct 
travel time = is) Obtained and from Equation (5), the reduced separation 
x Thus, the effects of the first layer are eliminated for each point and 
the problem reduced to a single layer caSe. 

The effect of dip is removed in a way Similar to that 


used for the first layer. Equation (17) becomes 





where 


and the subscript 'a' indicates the assumed values of the trial solution 


of Vo and 8. The reduced times corresponding to the flat layer case T, 
r 


30 





are obtained from 


ee ary (36) 


The reduced values are then used to compute the co- 


efficients of a straight line equation of the form 


2 
= 1p 
"2 Re k 


2 
Vy io 7) 


where Me should be zero because wW in Equation (36) eliminates T. 








ae 
and 
a ] 
2 
and Vo can be readily computed. The dip g,, is then corrected by 
a 
Vo 
tan 8, = i tan 8, : (3 9) 
a. a 


The solutions for Vo and 8, then replace the original 
trial solutions and the process is iterated a second time. A final 


solution is obtained for V,, 5, and H, and the computations move on 


Z Z 


to the next layer. 


Thus, in the first layer, knowing the sounding speed 


V. and the dip 8 


j , and having the record of reflected travel times T 


L t 


verses direct travel times Di, the value of the water depth Hy and 


ou 





horizontal sound speed V, were computed. In the second layer, the 


h 


problem was reduced to finding the interval sound speed Vo and the 
thickness H, by eliminating the effects of the first layer and then re- 
ducing the effect of dip uSing a trial solution. In each successive layer, 
the problem is first reduced to a single layer case by eliminating the 
effects of all upper layers using the computed values of thicknesses 
and interval sound speeds of those layers, and then the effect of dip is 
subtracted using an iterative process and a trial solution. 

This method is both fast and accurate. The time to 
reduce the data from a typical station is usually less than 20 seconds 
on the IBM 360/67 digital computer and the theoretical accuracy obtained 
is better than 1 part in 10,000 if the proper dips are entered (Le Pichon, 
Ewing, and Houtz, 1968). The primary advantage of this method is that 


there is no requirement that the layers be parallel, but rather, they may 


be dipping with respect to one another. 


OF Pane Oil U PER PROGRAM 

The HONDO I program used to reduce the wide-angle reflections 
is listed in Appendix C. It was adapted from a similar program used at. 
Lamont-Doherty Geological Observatory (LGO) of Columbia University 
(Houtz, personal communication) and follows the mathematical technique 
just presented of solving for thickness and interval sound speed of 
successive sub-bottom layers. Certain modifications were necessary 


to make the program compatible with the wide-angle reflection system 


2 





developed and for use in conjunction with the IBM 360/67 digital 
computer at the Naval Postgraduate School (NPS) using FORTRAN IV 
language. A thorough description of the routine followed by the program 
is presented in Appendix B. 

Sample data were included with the LGO program and modified 
accordingly (Appendix D) to be used with the HONDO I program. The 
corresponding results (Appendix E) were compared with those obtained 
using the IBM 1130 digital computer at LGO and the original form of the 
program, SLOW I. No significant difference was observed between the 
two solutions. The results from both programs are given in Table I. 

i. Data Reduction 

Instructions for entering the data on IBM cards are included 
on the first page of the HONDO I program list. A more delicate part of 
the data reduction consists of selecting and tracing the major reflectors 
from the precision seismic graphic record made during the field profiling. 
Here special care must be taken to avoid multiples that can occur and to 
use only those traces representing the events from true lithological re- 
flectors. The difference between true reflectors and multiples is shown 
in Figure 5. It was assumed that the interval sound speed increased 
with depth. The apparent reflection times from multiples at a given 
separation are greater than those of the corresponding true reflector 
because the multiple travels longer in the lower speed layer than the 
true reflector. Note also, the T/X trace of the multiple is nearly parallel 


to the true reflection from the same reflector. Results that are not 
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Comparison of Solutions of SLOW I 
and HONDO I Computer Programs 


Layer 


FIRST LAYER: 


Reflection Time 
(seconds) 
Sound Speed 
(km/sec) 
Propagated Error 
(km/sec) 
Depth 
(meters) 


SECOND LAYER: 


Reflection Time 
(seconds) 
Sound Speed 
(km/sec) 
Propagated Error 
(km/sec) 
Thickness 
(meters) 


PRD LAYER: 


Reflection Time 
(seconds) 
Sound Speed 
(km/sec) 
Propagated Error 
(km/sec) 
Thickness 
(meters) 


Table I 
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Figure 5. Multiple reflections 
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consistent with the local geology present another indication that multiples 
were present and traced. After the true reflectors are determined, the 
wide-angle reflection record is digitized in accordance with the pro- 
cedure outlined in Appendix A. 
Zs. Trial Solutions 

As discussed previously, the trial solutions are obtained 
from normal profiling records. The trial values of interval sound speed 
and dip for each Successive layer are entered on IBM cards in accordance 
with the instructions listed in Appendix C. If the interfaces are hori- 
zontal, there are no corrections for dip in the computations and the values 
for interval sound speed are unimportant. If the interfaces are not hori- 
zontal, inaccurate estimates of the sound speeds lead to only small 
errors in the first iteration and become negligible by the second iteration. 
However, dips entered which are inconsistent with the trial interval sound 
speeds may result in large errors. The error caused by entering such 
inconsistent values was investigated and the results given in Table II. 
The true solutions were those values computed using the data as given 
in Appendix D. The values computed in Case A were obtained assuming 
a trial sound speed 0.1 km/sec greater than the given data for both ne 
second and third layers. For Case B, the trial sound speeds were un- 
changed, but the dips of the Ac and third layers were entered as 


O 


|) It Can be seen from Table II that the error introduced by incorrect 


sound speeds is small, but entering dips inconsistent with the trial sound 
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speeds causes significant error. Furthermore, the computer program 
will not detect such errors but will compute values consistent with the 


dips entered. 


Table II 


Compeamson of Results of 
Various Trial Solutions 


Computer Sound Speed (km/sec) Propagated Error (+km/sec) 


Second layer Third layer second layer Third layer 
True Solution: i 7 ool Mer ors 0.06 0.10 
Case A: Le 7oot WG 9 0.06 0.10 
Case B: ie 7 Ol lenovo Or. O03 


The total time to compute the results (Appendix E) from the sample 
data (Appendix D) was 16.36 seconds on the IBM 360/67 digital computer. 
The accuracy of the results is considered to be better than 1 part in 10,000 
if the correct dips are entered. The propagated error gives no indication 
that the correct dips were entered or that the computed values represent 
the true interval sound speed. Only by subjectively analyzing the results 


can gross errors be realized. 
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IV. INSTRUMENTATION 


A. ibe ANGE kr LECTION SYSTEM 
The system for recording wide-angle reflection data is shown 
schematically in Figure 6. The profile recorder, timer switch, sound 
source, amplifiers and filters are those components normally carried 
on board oceanographic survey vessels for collecting normal-reflection 
seismic profiles. For collecting wide-angle seismic reflections, a 
sonobuoy system was developed to be used in conjunction with the on- 
board components. The sonobuoy system included a radio sonobuoy, an 
antenna, and a radio receiver. The system was developed with emphasis 
on simplicity and low cost; however, the two were not always Synonymous. 
eg The Sonobuoy 

The sonobuoy, launched while the ship was underway, per- 
formed as a self-contained station for detecting the sound pulses reaching 
that point in the ocean by various paths as shown in Figure 1. The 
pressure pulses were received by the hydrophone suspended below the 
sonobuoy by a thin wire. The pulses were converted to electrical signals 
and transmitted to the ship by the VHF radio transmitter contained in the 
watertight compartment of the sonobuoy. The sonobuoy transmitted on a 
pre-set frequency indicated by the channel number painted on the case. 
Table III lists the 31 sonobuoy channels available and the corresponding 
frequencies. The transmitter and acoustical amplifiers are powered by 


the sea-water-activated batteries located in the base of the housing 
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and become active within seconds after the sonobuoy contacts the water. 
The sonobuoys used were of the AN/SSQ-23A, AN/SSQ-41, and 
AN/SSQ-57 types. Figure 7 is a drawing of an AN/SSQ-57 sonobuoy 

as it would appear in the water showing the antenna and hydrophone 
package. The characteristics of the various types of sonobuoys used 


are listed in Table IV. 


Table III 


Sonobuoy Frequencies 


Channel Frequency Channel Frequency 
MHz MHz 
} hOZ.20 17 O7702 5 
2 es 2.010 18 632370 
3 Fos. 19 19 164.125 
4 iey4es Sil0) 20 P6407 5 
5 i Sre2 5 ih OOO 2a 
6 166.00 ae. 6G .0.75 
7 1667 5 ZO NG 7.125 
8 NOU 24 Ov O72 
9 16825 295 163.625 
io 169.00 26* 169.375 
el 169.75 Ze AO 125 
12* LO eo) Zoe 12067 5 
13 lal. 25 Z9 lw b2 5 
14 172.00 5 Or 7235 7S 
1S NET Zsoh ones pisos lo 
te Piso 


* crystals available for radio receiver 
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Figure 7. AN/SSQ-57 Sonobuoy 
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Table IV 


Sonobuoy Characteristics 


AN/SSQ-23A AN/SSQ-41 AN/SSQ-57 
Weight, 8.11 9.02 9.02 
Kg 
Channels 1-16 1 - 3l 1 - 3l 
Time Selection, l il one 4 Orc 
Hours 
Depth of Hydrophone, S02. Sz Or ro.) 76.9 
Meters 
Transmitter Power, 72 l | 
Watts 


Each sonobuoy has a water-soluble scuttling plug located 
above the center of the housing. When the plug dissolves, water enters 
the watertight compartment and the buoy sinks. The floating period of 
the sonobuoy is a function of water temperature, with the plug dissolving 
more rapidly in warmer water. In no case, however, is the floating 
period less than 8 hr nor more than 20 hr in the case of the AN/SSQ-57 
type (Naval Air Systems Command, 1967). 

The plug was sealed with waterproof tape when a study was 
conducted investigating modifications to be made to the sonobuoy for 
recovering and re-using it on subsequent stations. The study revealed, 
however, that the rubber gaskets used to seal the water tight compartment 


were not intended for long use and that small amounts of sea water entered 


42 





the compartment after only 3 or 4 hr of immersion. This was learned after 

a considerable amount of sea water was found inside one particular 

sonobuoy after it was in the water for about 12 hr. On another occasion, 

a sonobuoy operated successfully for 3 hr in the water before it was 

recovered. The salt water batteries were removed immediately to prevent 

any electrolysis or heating damage and the buoy was washed ina mild 

solution of soap and water. Later, however, it was found to be inoperable 

due to corrosion of the electronic components inside the watertight case. 

In each case, the scuttling plugs were sealed. Considering this, and 

also considering the problem of trying to recover 18 to 90 m of hydrophone 

wire, which had a strong tendency to kink, it was decided that the cost 

of using one sonobuoy per station was far less than the complications 

which would incur and the time involved to recover and re-use a sonobuoy. 
The sonobuoys were launched by hand from the fantail of a 

ship. A screwdriver was used to bend the ends of the retaining ring on 

Dre Dottomeartme Sonobuoy and release the bottom plate (Fig. 8). In 

some cases, the small screw used to secure the ring was removed and 

the entire retaining ring removed. Once the ring was removed, the 

bottom plate and hydrophone package were free to drop out of the sonobuoy, 

but were held in place until the sonobuoy was launched. The bottom 

plate serves as a weight to pull an aluminum protective housing away 

from the hydrophone package once the sonobuoy is in the water. For this 

reason, the clips (Fig. 9), which secure the plate to the housing, are not 


femoved . 
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Figure 8. Releasing the bottom plate of the 
AN/SSQ-57 sonobuoy 
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Figure 9. Bottom plate retaining clips 
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The next step was to release a second retainer ring 
located around and near the top of the sonobuoy. Considerable caution 
was warranted here as this ring allowed the rotochute assembly to fall 
free and pull a small white lanyard releasing the spring-loaded transmitting 
antenna. The antenna erects with sufficient force to cause injury. Once 
it was safe to do so, the ere ring was released by lifting a small 
U-shaped pin (Fig. 10) and the rotochute assembly gently removed. A 
slight tug given on the lanyard tripped the antenna retainer clip allowing 
the antenna to extend to its full operating position shown in Figure 7. 

The sonobuoy was then ready to be launched. 
Zs Receiver System 

A radio receiver and antenna system was used to receive 
and amplify the VHF signals transmitted by the sonobuoy. The radio 
used for this was a Japanese-made Realistic Patrolman pro-2 solid state 
VHF receiver. The decision to use this radio was based mainly on its 
availability at the time the study was started and the fact that its frequency 
band covered that of the entire sonobuoy range. Also of importance, and 
in keeping with the idea of simplicity and economy, the pro-2 was small, 
easy to carry, Simple to operate, and relatively inexpensive ($99.00). 

The pro-2 radio receiver is a self-contained unit operated 
Srenermal tlO-V AC. An internal speaker was used to monitor the signals. 
The signal was tuned by manual control for maximum volume. Crystal 
operation was available, but did not appear to increase the quality of 
the signal. The crystal did, however, provide a means for tuning the 
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receiver without a signal present. With the proper crystal in place, the 
radio was simply tuned for maximum background noise. Crystals for the 
seven sonobuoy channels indicated in Table III were purchased from 
Radio Shack Warehouse, 2615 West 7th Street, Fort Worth, Texas. 

The pro-2 receiver had two locations where the signal was 
available for recording. The signal taken from the speaker jack on the 
front panel of the cabinet was found generally to be highly distorted by 
the audio amplifier section of the receiver. The tape-out jack on the 
back of the cabinet provided a signal of rather low voltage, but when 
amplified it was found to be relatively undistorted. The results were 
further verified by using an oscilloscope to observe the wave form of 
the signal from each jack. A pure tone was transmitted froma signal 
generator and received by the radio. The signal at the tape-out jack 
was clear and undistorted while that of the speaker jack was either 
distorted at high volume settings or was lost in the amplifier noise at 
low settings. 

An ll-element Winegard Yagi antenna ee 11) was used 
with the pro-2 radio receiver for receiving the signal transmitted by the 
sonobuoy. The Yagi was purchased when the ground plane dipole antenna 
(salvaged from a used sonobuoy) failed to provide sufficient power gain 
for receiving the signal from the sonobuoy except at very small separation 


distances. 
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Figure 11. Yagi ll-element antenna 
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The antenna patterns for each antenna were recorded over a 
range of frequencies covering the sonobuoy range (Fig. 12). It can be 
seen from this graph that the Yagi was quite frequency-dependent and 
at lower frequencies did not produce better results than were obtained 
using the dipole antenna. Figure 13 shows the antenna pattern for the 
Yagi at various frequencies. The dipole antenna pattern was omni- 
directional over the entire sonobuoy frequency range. It was decided, 
therefore, that only higher frequency sonobuoys would be used and 
special attention would be given to insure that the Yagi antenna was 
pointed in the direction of the sonobuoy. 

A 300/50-ohm impedance matching transformer was attached 
at the antenna to couple the 300-ohm antenna to the 50-ohm shielded 
cable. All external connections were covered with alternate layers of 
tape and electrical coating to prevent salt spray from damaging the system. 
A stainless steel pipe, 7-m long and 3.75 cm in diameter, was used for 
an antenna support. 

The entire wide-angle reflection system weighed less than 
10 kg per sonobuoy and was quite transportable; the whole system can 
be carried in the backseat of an automobile. It took only a few minutes 
Pomsetup tie equipment aoeard Ship. The system was simple to use, 


requiring little more than to launch a sonobuoy and to tune the receiver. 
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Figure 13. Antenna patterns for Yagi antenna at 
various frequencies (relative to 175 MHz) 
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Vee CONCLUSIONS 


It was found from this study that the Realistic pro-2 radio receiver, 
when used with a Winegard Yagi antenna, could adequately receive the 
VHF signal transmitted by a radio sonobuoy and at sufficient distance to 
conduct wide-angle seismic reflection profiling. It was further found 
that the signal received with this system was distorted by the audio 
amplifier section of the radio. An undistorted signal was available at the 
tape-out connection on the back panel of the pro-2 receiver cabinet, by- 
passing the audio amplifier section. 

It was concluded after becoming familiar with three types of radio 
sonobuoys that, when used normally, they were quite simple to operate; 
but any modifications that were possible required a considerable amount 
of preparation in order to use the sonobuoy for more than one station. 

Magnetic tape recordings have been made of the signal received 
using the sonobuoy system and the sample set of data from LGO demon- 
strated that the HONDO I computer program adapted for use with the 
system was able to reduce wide-angle reflections and obtain the interval 
sound speeds and corresponding thicknesses of sub-bottom layers. 

It is concluded, therefore, that the system developed can be used 
to collect wide-angle seismic reflection data and that the HONDO I 


program can be used to reduce the data. 
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ee SIND TX, AK 


DiC iZzinieetHh WIDE-ANGLE REFLECTION RECORD 


If TS is marked on the record (the instant when the sonobuoy and 
sound source are together), this point becomes the origin. Otherwise, 
the direct travel time line, called the direct line, is extended to intercept 
the X-axis, thus defining the origin. A plastic overlay is placed over 
the record so as to include the origin and all direct and reflected traces 
within the millimeter-size grid etched on it. 

Figure 14 shows a sample record with a direct line and one reflected 
trace. The origin of the grid does not coincide with the origin of the 
record, but rather the entire data is included under the grid. Below the 
record are sample values read from the record, including the two points 
A and B used to set the time scale. Points A and B are selected anywhere 
along the T-coordinate and are a measure of the distance in millimeters 
from the origin of the grid in the T-direction only. The time DDD rep- 
resents the known time interval between the two points. 

The values from the direct line are recorded in eight-digit integer 
form. For example, a point located at 127.9 mm from the grid origin on 
the X-scale and 148.9 mm on the T-scale is recorded as 12791489. The 
X-scale is the vertical scale and represents the first four digits, while 
the T-scale is horizontal and represents the last four. From two to six 
points can be used to represent the direct line. The reflection traces 
are digitized similarly, except that as many as 90 points can be used. 


The program can handle up to 10 reflecting horizons. 
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Figure 14. Digitizing the data 
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At times, the sound pulse is detected at the sonobuoy after the 
recorder completes a cycle. This occurs ee the separation between 
the sonobuoy and the sound source oe large enough so that the 
travel time for a particular sound wave is greater than the sweep rate 
of the recorder. An example of this is shown in Figure 14 where the 
trace goes off the bottom of the record and continues on the top during 
the next cycle. In this case, the shift is taken care of within the 


program. 
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APPENDIX B 


ROUTINE FOLLOWED BY HONDO I COMPUTER PROGRAM 


The computations begin in the MAIN routine by reading from the 
first card in the data deck the station identification, ADEN] and ADEN2, 
and the number of layers or reflections observed, NREF, for the particular 
Station. Subroutine TRIALS is then called to enter the valves for the trial 
solutions which include trial sound speeds V from the second card in the 
data deck and the trial dips W from the third card. Only NREF values are 
taken from each card and blank spaces are read as zero. AS an example, 
for the sample station DIANNE (Appendix D), the card representing the 
trial dips is blank; therefore, the trial dip for each of the three layers 
is zero. The technique of reading zeros from blank cards is used through- 
out the program to control the exits from particular subroutines. For this 
reason, it is extremely important that a blank card follow the cards con- 
rman tne digitized data pairs for the direct trace and for each reilected 
trace, and that an additional blank card follow the data from the last 
layer of the last station investigated. This additional blank card causes 
the program to exit the computer and completes the job. The blank card 
after each set of data points tells the computer that there are no more 
points for that particular trace. 

ive tounehec hemmetiercatesceck 16 also read by TRIALS. The card 
contains the value of the computed sounding speed VV for the water layer. 


This value will replace the value read in as a trial speed for this layer. 


of 





Thus, the first number in the array labeled V is the actual value used for 
the sounding speed VV of the first, or water, layer. Also, the dip read 
for the first layer is not, in fact, a trial dip; but rather it represents the 
actual value determined for the sea floor. For this reason, the values 
entered for the first layer of sounding speed VV and dip W must be those 
values computed from observations made at the time the data were col- 
lected. 

Values for V and W are transferred back to MAIN through COMMON. 

The program is returned to MAIN and the main loop for NREF layers 
begins by setting the integer L equal to | for the first layer. MAIN then 
calls both subroutines READS and INTERP in succession. 

For the first layer, READS enters the values from the fifth data card 
representing the constants which appear in Figure 14. They include two 
arbitrary points on the time scale A and B, the time in seconds represented 
by the distance between them (DDD), the seconds required for each sweep 
of the recorder (SPS), and finally, DELAY, which is a provision to include 
any time Atlee which might occur between the shot and the start of the 
recorder sweep. READS uses A, B, and DDD to scale the time axis in 
seconds. SPS is used to add the appropriate number of seconds to each 
point of the trace if the trace is continued at the top of the record during 
a scale shift. 

Once the time axis is scaled, READS then enters the digitized data 
points from the direct travel time trace. These points are converted to DX 


in millimeters and DY in seconds from the origin. 
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Next, READS enters the digitized points from the first reflected 
time trace and converts these points to RX and RY (in millimeters and 
seconds, respectively, from the origin). 

Subroutine COUNT is called from subroutine READS to count the 
number of non-zero data points found on each card. The total number of 
points for any reflected time trace is summed on II and carried via 
COMMON to other locations in the program. ND is the number of points 
used to compute the direct time line. 

All values entered in READS are transferred to MAIN and then to 
INTERP through COMMON. 

INTERP interprets RX in millimeters for DD in seconds of direct 
travel time corresponding to each reflected travel time RY. This is 
accomplished by computing the slope of the direct time line, DELXY, 
and multiplying it by the millimeters of separation between RX and the 
origin DX1 as shown in Figure 15.RX is then reset equal to DD at each 
point on the reflected trace. Thus, RX becomes the travel time corresponding 
to each penieeted travel time RY. RX and RY are transferred to MAIN 
through COMMON as U and S respectively. In MAIN, the direct travel 
time corresponding to the first reflection time recorded is compared with 
the minimum allowable separation for considering a fourth-order polynomial 
fit to the data, DMIN. DMIN is equal to 0.3 sec. If the direct time is 
less than DMIN, subroutine LSFIT is called and the data are used to 
compute a fourth-order polynomial using a least-squares technique. If 


the data are not obtained at or near vertical incidence and the minimum 


Be. 





DX | RX 





DD = DYI1+ (RX — DX1) * DELXY 


Figure 15. Relating direct travel time to 
corresponding reflected travel time 
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reflection time is greater than DMIN, the data are squared and a linear 
least-square polynomial is computed in subroutine LSFIT. In either case, 
the data are extrapolated to Bo using the resulting polynomial equation,,. 
the minimum reflection time is obtained, and the depth is computed using 
the sounding speed. 

Subroutine LSFIT also computes the theoretical reflection time 
corresponding to the measured direct time by solving for the equation of 
a line fitted to the th? data.. 

Back in the MAIN routine, the corrective term CON to be applied 
to the travel time is computed and subtracted from the total time for 
each point. Thus, the reduced times R and D are squared and LSFIT is 
called to fit the values to a linear least-squares line and compute the 
horizontal sound speed from the coefficient X(2). 

This completes the computations for the first layer, having obtained 
the depth and horizontal sound speed. These results are printed and the 
program returns setting L equal to 2. 

In successive layers, the program begins by calling the subroutine 
READS. The digitized data for the particular layer is read and converted 
to millimeters of separation and seconds of travel time from the origin. 
Subroutine INTERP is called, interpreting millimeters of separation for 
seconds of direct travel time as was done in the case of the first layer. 
Now, however, the horizontal sound speed is known and the separation 


can be converted to distance in kilometers by an equation equivalent to 
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Equation (5). The values of reflected times and corresponding values of 
horizontal separation are returned to MAIN through COMMON as S and U 
respectively. 

In MAIN, the horizontal separation for the first point is compared 
with the minimum allowable separation for considering a fourth-order 
polynomial, DMIN. In this case, DMIN is 0.45 km, obtained from the 
product of an assumed horizontal sound speed of 1.5 km/sec anda 
minimum allowable separation, time of 0.3 sec. If the first point is less 
than DMIN, the data are used in the subroutine LSFIT to compute the co- 
efficients of a fourth-order least-squares polynomial. If the first point 
is greater than DMIN, the values are squared and used to compute a 
linear lease-squares equation, again in subroutine LSFIT. As was done 
for the first layer, the data are extrapolated to obtain the minimum re- 
flection time for the particular layer being considered. 

The derivative of the fourth-order polynomial with respect to the 
horizontal separation is used to compute the angle of emergence 
ZM (8 + 6 A Fig. 4) using an equation equivalent to Equation (33). If 
the minimum reflection time was previously computed using a linear fit 
to the data, subroutine LSFIT is called a second time to compute the 
coefficients of the fourth-order polynomial fit. ZM is a function of the 
slope (dT/dX) of the fourth-order curve. The curve is not restrained 
beyond the first and last points and its slope may become erratic there. 
For this reason, the angle of emergence is not computed for the first 
and last points. 
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Subroutine TUPPER is called and the thickness of the layer being 
considered is computed using the trial solution of interval sound speed 
and dip. The technique follows that given in the previous section for 
the first layer, using an equation equivalent to Equation (29). 

Subroutine TUPPER uses the emergence angles to compute travel 
times TT in the layers above the layer for which a solution is sought. 

TT is computed using the previously-solved values of interval sound 
speed and dip for the upper layers. The corresponding wales of hori- 
zontal separation are found ina similar way. The travel times and 
separations are transferred to MAIN as TR and DR respectively. 

MAIN reduces the measured values of reflection time S and 
horizontal separation U to obtain D, the reduced value of horizontal 
separation, and R, the reduced reflected travel time. Dis also adjusted 
for the sum of the dips of the upper layers by dividing by C, the cosine 
of the sum. 

The reduced D and R corresponds to a single layer case as shown 
imenigure 3 . 

The travel time for each point is finally reduced to that of a flat- 
layer case by subtracting equations equivalent to Equations (10) and (15). 
The value R in the computer program is now the reduced time squared for 
the flat-layer case. The value D is also squared and the subroutine 
LSFIT is called to fit a line to the reduced values. The second coef- 
ficient X(2) of the linear equation obtained from the subroutine LSFIT is 
then used to compute the interval sound speed. 
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The dip of the layer W is corrected using the relation 


tan W = - tan W_ , 
V a 


a 
where the subscript ‘a' indicates the values for the assumed trial 
solution. 

This solution then replaces the original trial solution and the 
computations go back to computing a new thickness and start the second 
iteration. Upon completion of the second iteration, the solutions for 
the layer are printed and the routine starts a new layer. 

The results are tabulated and printed for each layer in succession. 
A sample list of results for the theoretical station DIANNE are included 
in Appendix E. 

The first line gives the station name, layer number, and number 
of points used by the programs to compute the results. The number of 
points considered may be less than are digitized due to rejection of 
points occurring before L. or those resulting in negative reduced times. 
The first three columns of figures below the heading represent squared 
values of reduced separation and reduced reflected travel times cor- 
responding to a flat-layer case. For the first layer, the square of the 
separation listed under the column labeled X2, are in seconds squared, 
whereas, for each of the following layers, they are in kilometers squared. 

The columns labeled T2 MEASURED and T2 COMPUTED are the 


reflection times squared corresponding to the separation in the first 
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column. The measured values are those values read from the reflection 
curves and reduced to the flat-layer case. The computed values are 
from a least-squares curve fitted to the reduced data. The fourth column 
is the difference between the measured and computed values of the 
reduced reflection times. 


Z 
The variance o of the data is computed from 


where N is the number of points, T the observed times, and T. the 
computed times. The standard deviation is the positive square root of 
the variance. 

If many reflection profiles were to be taken at the same location, 
the variance of the interval sound speed about the mean can be calculated. 
However, with only a single set of data, the variance cannot be calculated 
in this way. Instead, the techniques used to compute the variance, or 
propagated error, for each interval sound speed is based on the variance 


of the coefficients for each point compared to those coefficients obtained 


with the least-Squares fit. 
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APPENDIX C. HONDO I COMPUTER PROGRAM 
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SUBROUTINE READS 
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APPENDIX E. RESULTS 
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